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A general technique is developed for calculating func- 
tional determinants of second-order differential operators with 
Dirichlet, periodic, and antiperiodic boundary conditions. As 
an example, we give simple formulas for a harmonic oscillator 
with an arbitrary time-dependent frequency. Here our result 
is a generalization of Gel'fand-Yaglom's famous formula which 
was restricted to Dirichlet boundary conditions. Apart from 
the generalization, our derivation is more transparent than 
theirs, the determinants requiring only knowledge of the clas- 
sical trajectories. Special properties of operators with a zero 
mode are exhibited. Our technique does not require the cal- 
culation of the spectrum and is as simple as Wronski's method 
for Green functions. 

1. Evaluation of Gaussian path integrals is necessary in 
many physical problems. In particular, it appears in all 
semi-classical calculations of fluctuating systems. Typ- 
ically, it leads to the problem of calculating the func- 
tional determinant of a second-order differential opera- 
tor pi. For Dirichlet boundary conditions, a first gen- 
eral solution of this problem was given by Gel'fand and 
Yaglom [pi based on the lattice approximation to path 
integrals in the continuum limit. Their result was ex- 
pressed in terms of a simple differential equation for the 
functional determinant. In subsequent work P]- M, the 
formalism was generalized to a variety of differential op- 
erators and boundary conditions based on the concept of 
zeta-function regularization 0|. Unfortunately, Gel'fand- 
Yaglom's method becomes rather complicated for pe- 
riodic and antiperiodic boundary condition relevant in 
quantum statistic (see Section 2.12 in |l|). In the periodic 
case there is, moreover, a zero mode causing additional 
complications. 

In this paper we present a systematic method for find- 
ing functional determinants of linear differential opera- 
tors which is based on Wronski's simple construction of 
Green functions. Our method is simpler than those used 
in the previous approaches, since the determinants are 
expressed entirely in terms of a classical trajectory. Fur- 
thermore, for fluctuation operator with a zero mode, a 
case frequently encountered in semiclassical calculations, 
the special treatment of this mode becomes transparent. 

2. The typical fluctuation action arising in semiclassical 
approximations has a quadratic Lagrangian of the form 



[i^ -n'^{t)x^ 



(1) 
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Physically, this Lagrangian describes a harmonic oscil- 
lator with a time-dependent frequency Q{t). The path 
integral for such a system was studied in several papers 
Q~ |ll|. For such an oscillator, both the quantum me- 
chanical propagator and the thermal partition function 
contain a phase factor exp[i5'ci(x)] and are multiplied by 
a pre-exponential fluctuation factor proportional to 



F{h,ta) 



/BetKi 
V DetK 



-1/2 



(2) 



where Ki = -d^ - n^{t) = Kq - n^{t) is the kernel of 
the second variation of action S{x) along the classical 
path Xci{t). The linear operator Ki acts on the space of 
twice differentiable functions y{t) = Sx{t) on an interval 
t S [ta,ti,] with appropriate boundary conditions. These 
arc Dirichlet boundary conditions y{ta) — y{tb) = 
in the quantum-mechanical case, and periodic (antiperi- 
odic) y{tb) — ±y{ta),y(tb) = rkyita) in the quantum sta- 
tistical case. In these two cases the operator K may be 



chosen as Kq or Kq 



■^0 5 



respectively, where ujq is a 



time-independent oscillator frequency. The ratio of de- 
terminants (0) arises naturally from the normalization of 
the path integral and is well-deflned [Q . Furthermore, for 
such an operator we may assume the Fredholm property 



DetJCi 
Det^ 



T>cik-^Ki 



(3) 



thus neglecting multiplicative anomalies |12| . Since the 
operator K^^Ki is of the form I + B, with B an operator 
of the trace class, it has a well-deflned determinant even 
without any regularization. 

To calculate F(tb,ta), we introduce a one-parameter 
family of operators Kg depending linearly on the param- 
eter g : Kg — Kq - gVi^ (t) , -C 5 ^ 1 . The above prop- 
erty (H) allows us to make use of the well-known formula 
logDetii'"^^^ = Tr\ogK~^Kn to relate the ^-derivative 
of the logarithm of the ratio (0) to the trace of the Green 
function of the operator Kg as follows 

dg \ogV)eik-^Kg = -Tr n'^{t)Gg{t, t'), (4) 

the Green function being deflned by 

Gg{t,t') = [-d?-gn\t)]-'Sit^t'). (5) 

Formula (0) is valid provided we regularize the trace 
on the right-hand side, if it diverges, via zeta-functions 



C(s) = Si'^""' where the sum runs over all eigenval- 
ues. It is convergent for sufficinetly large s and defined 
for smaller s by analytic continuation (see 0). Then, 
for each member of the g-family, DetKi = exp[— ^'(0)]. 
Another proof of fl) can be found in . 

By integrating (^, we obtain for the ratio of functional 
determinants (0): 



BctK-^Kg^Ce^pl- dg' dtn^{t)Gg'{t,t) 



(6) 

where C = DetK~^Ko is a ^-independent constant. This 
is our basic formula to be supplemented by an appropri- 
ate boundary condition to Eq. (0) for the Green function 
as we shall now discuss in detail. 

3. A general solution of Eq. (^) is given by advanced 
or retarded Green functions as follows 



G;it,t')=G+it',t) = eu'-f9it,t'), 



(7) 



where Qtt' = 0(t — i') is Heaviside's function and fg{t, t') 
is a combination 



fait^t') 



1 



[VgmAt') - Ut)Vg{t')] 



(8) 



of two linearly independent solutions rig{t) and ^g{t) of 
the homogeneous equation 



[-d^--gn^{t)]hg{t)^o. 



(9) 



The constant Wg is the time independent Wronski de- 
terminant Wg = rjg^g " f]gS,g- Thc solutlou (|7|) is not 
unique since it leaves room for an additional general so- 
lution of the homogeneous equation (y) with an arbi- 
trary coefficients. This freedom is removed by appro- 
priate boundary conditions. Consider first the quantum 
mechanical case which requires Dirichlet boundary condi- 
tions y{th) — y{ta) — for the eigenfunctions y{t) of Ki, 
implying for the Green function the boundary conditions 



Ggita,t')=0, t<t'. 

Gg{t',tb)^0, t' <t. 



(10) 



The operator K in the ratio (||) is equal to Kq, and the 
constant C in Eq. (||) is unity. After imposing ([lO|), the 
Green function is uniquely given by Wronski's formula: 



Gg{t,t') 



where 



Qtt'fgit', ta)fg{tb, t) + Qt'tfg{t, ta)fg{tb, t') 



fgita,tb) 



fg{ta,tb) 



DetAg 

w„ 



^0, 



(11) 



(12) 



with A being a constant (2 x 2)-matrix 

Vb 6 



(13) 



formed from the solutions r]g{t) and ^g{t) at arbitrary 
5 7^ 1. Note that these solutions are restricted only the 
condition (p^). The result is unique and well-defined, 
assuming the absence of a zero mode ^(t) of the operator 
Ki with Dirichlet boundary conditions £,a = ^b = 0. Such 
a mode would cause problems since according to ( [lO| ) , the 
Wronski determinant W would vanish at the initial point, 
and thus for all t. 

Excluding zero modes, we obtain from (pi): 



Trn^{t)Gg{t,t')^ 



/ f ' dtn^t)fg{t,ta)fg{tb,t). 

Jg(ta,tb)Jt^ 

(14) 



To perform the time integral on the right hand side, we 
make use of the identity 

nHt)mvit) = dt[Ut)d9Ut) ~ Vgit)dgUt)]. (15) 

This follows from Eq. (||) for r]{t,g), and an analogous 
equation for ^g (t) , after multiplying the first by ^^ (t) and 
the second by T]g{t), and taking their difference. In the 
limit g — + q, we obtain (|l5|) from the linear term in g — g. 
Inserting ( p5| ) into (|lj), we see that 



Trn^it)Gg{t,t')^~dglOg(^^^ 

Substituting (Rq) into (H), we find 

DctAg /DetAo 
'W^/ Wo ' 



BctKa^Kg ^ 



(16) 



(17) 



where DetAo/ Wq = tb — ta- Finally, setting to g = 1 in 
( [l7| ) gives the required ratio of the functional determi- 
nants 



BetK,'-K,^'^^/(tb-ta). 



(18) 



In a time-sliced quantum mechanical path integral, the 
determinant of Kq is finite and has the value M 



BetKo = tb-ta, 



so that we obtain 



DctATi = 



[Vl{ta)^l{tb) - VlitbMtq)] 

Wi 



(19) 



(20) 



which coincides with Gel'fand-Yaglom's formula (see Sec- 
tion 2.7 in g). 

For a harmonic oscillator with a time-dependent fre- 
quency f2(i) it is convenient to relate the set of two in- 
dependent solutions r]g(t) and ^git) of Eq. (||) at g = 1, 



for which we omit the subscripts g, to the classical path 
a^ci(^) = Xa^{t)+xi,r]{t) satisfying the endpoint conditions 
Xci{ta) = Xa and Xci(ifc) — xi,. Since this construction 
satisfies ?7q = 6 = 0, Ty^ = ^a = 1 and W = ^b = -7?a, 
the explicit solution being 



m 

rjit) 



dxcijt) _ p{t)pbS[nuJo{qb - q) 
dxa PaPb sin ujoiqb- qa)' 

dxcijt) _ p{t)paSmuJo{q -qg) 
dxb PaPb sin ujo{qb - qa) ' 



(21) 



They are parametrized by two functions q{t) and p(t) 
satisfying the constraint 



ujaqp = 1, 



(22) 



where wq is an arbitrary constant frequency. The func- 
tion p{t) satisfies the Ermakov-Pinney equation |13| ] 



p + n\t)p~p-^ = 



(23) 



Inserting (|2|) into (jlj), we obtain for the harmonic os- 
cillator with a time-dependent frequency r2(i) the ratio 
of functional determinants 



V)ciK^^Ki = 



PaPb sin uJo{qb - qg) 
{h - to) 



(24) 



where subscripts a and h indicate avaluation at t — tg 
and t = tb, respectively. We check this representation 
by expressing the right-hand side in terms of the clas- 
sical action S'ci(a;). With the same normalization as in 
(|l9|), this yields the well-known one-dimensional Van- 
Vleck formula 



DetKi 



~M[d'^Sci{xg,Xb)/dxadxb] 



(25) 



To end this section we note that the ratio (Esl) can 
easily be extended to the stochastic case where the final 
position of the trajectory xjt) remains unspecified. To 
this end we consider Eqs. dl4) and (Hq) with a variable 
upper time t' > t > tg. Then the eigenvalues of the op- 
erator K^^Ki become functions of t' with a phase factor 
produced by each passage through a focal point. 

4. Consider now periodic (antiperiodic) boundary con- 
ditions y{tb) = ±y(ia), y{tb) = ±y{tg) for the eigenfunc- 
tions y{t) of the operator Ki and the for Green function 

Ga{t,t'): 



G^tb,t')^±G^tg,t'), 
G°(4,t') = ±G°(ia,0, 



(26) 



where T — tb — tg is the period. In both cases, the fre- 
quency Q{t) and Dirac's (5- function in Eq. (||) are also as- 
sumed to be periodic (antiperiodic) with the same period. 
The general solution of Eq. (M) satisfying the boundary 



conditions (^6|) is constructed by adding to (0) an ex- 
pression of the same type as before, using the same ho- 
mogeneous solutions r]g{t) and ^g{t). The result has the 
form 



G;{t,t')=Gg{t,t') 

\fg{t,tg)±fg{tb,t)][fg{t',tg)±fg{tb,t')] 



A° -fgita^tb) 



(27) 



with the condition 






where A" are now the (2 x 2)-constant matrices 

iVbTVa) (6TCa) 



A" = 



{VbTVa) (CbT^a), 



(28) 



(29) 



evaluated at g ^ 1. In analogy to Eq. (|16|) we now find 
the formula 



Trn'it)Gl{t,t')^~dglogl^^ 



(30) 



Substituting this into (m and setting g — 1, we obtain the 
ratio of the functional determinants for periodic bound- 
ary conditions 

Deti^-i . K, = ^/4sin^ 5^V4_0 (3,^ 

Here Det K = Det(— 9j^ — ujq) is the fluctuation determi- 
nant of the harmonic oscillator, which in the same nor- 
malization as in (|l9) is equal to 



-n. i r> A ■ 2 ^0(4 - ta) 

Det K = A sm — ^ , 

2 

and thus the formula 

{rib - Va)iib - ia) - (6 - ia){rib - Va) 



BetKi = 



W 



(32) 



, (33) 



the right-hand side being evaluated a,t g — 1. For an- 
tiperiodic boundary conditions, the analogous expres- 
sions are 

DetX,- = ^Acos^ Mtb_tg^^ (34) 



p^^^ ^ iVb + ??a)(6 + Ca) - (6 + ^a){Vb + Va) ,^^. 

For a harmonic oscillator with a time-dependent fre- 
quency rj(i), we use again the representation ( pi| ) for ^(i) 
and ri{t) in terms of the functions p{t) and g(i)) which in 



addition to (|2^) and (g3|) have the following properties: 
the function p(t) is periodic and even 



p{t + T)=p{t), p{^t)^p{t) 



(36) 



so that pb — Pa, whereas the function q[t) satisfies 

q{t + T)=q{t)+qi,, qa = 0, (37) 

where T = {tb — ta). Inserting now the solutions (|2l| ) 
into IBWj and (M), we find the ratio of functional deter- 
minants for a harmonic oscillator with a time-dependent 
frequency fl(t) with periodic boundary conditions 



Betk-^Ki =4sin' 



2 ^^096 



/4sm — , 



and with antiperiodic boundary conditions 

DetK-^K, = 4cos2 ^^^Mcos^ ^. 
2 ' 2 



(38) 



(39) 



Note that only formula (|2J) for the Dirichlet boundary 
condition has been known in the literature (see Isl - f O | ) . 
The periodic and antiperiodic formulas (38) and (|39|) are 
new, although they have had predecssors on the lattice 
[ fi5[ . Moreover, our new derivation has the advantage 
of employing only Wronski's simple construction method 
for Green functions. The general expressions for the func- 
tional determinants (|20|), ( |33| ) and (|35| ) are form-invariant 
under an arbitrary changes (77, ^) —>■ {fj, ^) of the basic set 
77 (t) and ^(t) of two independent solutions of the homo- 
geneous equation (||). 

5. Contrary to the case of a harmonic oscillator with 
a time-dependent frequency 51 (t), consider now the situ- 
ation where the operator Ki has a zero mode. In this case 
we may assume the frequency fl{t) in Eq. (|l|) the special 
form n'^{t) = V"{xci{t))/M with a potential V{x), al- 
lowing reflecting the translation invariance of the theory 
with Lagrangian (m along the time axis flO]. Let f(t) 
be the corresponding eigenfunction satisfy the condition 
^a = as well as ^f, = 0. As mentioned above, the con- 
dition (|lj) is now violated, making Eq. (^l]) undefined, 
and it is impossible to construct two independent solu- 
tion ^(t) and ry(i) since their Wronski determinant would 
be equal to zero indentically W — ija^a — VaCa = due 
to the boundary conditions (10). Since the Wronski con- 
struction is not applicable we replace 



6=0, Ca = 

by the regularized conditions 

a=0, ^'a = e. 



(40) 



(41) 



These do not require a new calculation of the determinant 
(vw, and we find immediately 



6 



(42) 



in the limit £ ^ 0. We therefore remove the zero mode 
from the determinant using the standard method [[l7[ . 
The regularized determinant is defined by 



DetiCf = hm -, 



(43) 



where A"^ is the eigenvalue associated with the eigenfunc- 
tion ^(t). 



Kie = A^r 



(44) 



with the limits ^^ — > ^ A^ ^ for e ^ 0. To first order 
in e it follows from (E3) that 



dt^K^e~>^' dte{t) = \'{m- 



(45) 



Integrating the left-hand side by parts and taking into 
account the conditions (Eoh and (Ell) gives 



my 



(46) 



Finally, substituting (^ and ( f42| ) into (|4^) we obtain 
the functional determinant without zero mode 



T>eiK{ 



R_<i\^> 



iaib 



(47) 



For periodic (antiperiodic) boundary conditions y{ta) 
±y{ti,),y{ta) — ±y(th), the analogous formula is 



VaiVaL -Va^b) 



(48) 



In the periodic case, formula ( [48|) is useful for scmiclas- 
sical calculations of path integrals processing nontrivial 
classical solutions such as solitons or instantons [y. 
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